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We calculate the azimuthal dependence of the heavy-quark-initiated 0(a s ) contributions to the 
lepton-nucleon deep inelastic scattering (DIS). It is shown that, contrary to the photon-gluon fusion 
(GF) component, the photon-quark scattering (QS) mechanism is practically cos 2<^-independent. 
We investigate the possibility to discriminate experimentally between the GF and QS contributions 
using their strongly different azimuthal distributions. Our analysis shows that the GF and QS 
predictions for the azimuthal cos 2tp asymmetry are quantitatively well defined in the fixed flavor 
number scheme: they are stable, both parametrically and perturbatively. We conclude that mea- 
surements of the azimuthal distributions at large Bjorken x could directly probe the intrinsic charm 
content of the proton. As to the variable flavor number schemes, the charm densities of the recent 
I— | CTEQ and MRST sets of parton distributions have a dramatic impact on the cos 2ip asymmetry in 

, the whole region of x and, for this reason, can easily be measured. 
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I. INTRODUCTION 

The notion of the intrinsic charm (IC) content of the proton has been introduced over 25 years ago in Refs [US]. It 
was shown that, in the light-cone Fock space picture 0, it is natural to expect a five-quark state contribution to the 
proton wave function. The probability to find in a nucleon the five-quark component \uudcc) is of higher twist since 
it scales as 1 /to 2 where to is the c-quark mass Q ■ This component can be generated by gg — > cc fluctuations inside 
the proton where the gluons are coupled to different valence quarks. Since all of the quarks tend to travel coherently 
at same rapidity in the \uudcc) bound state, the heaviest constituents carry the largest momentum fraction. For this 
reason, one would expect that the intrinsic charm component to be dominate the c -quark production cross sections 
■ at sufficiently large Bjorken x. So, the original concept of the charm density in the proton [l|, || has nonperturbative 
Oh! nature and will be referred to in the present paper as nonperturbative IC. 

^ . A decade ago another point of view on the charm content of the proton has been proposed in the framework of the 
variable flavor number scheme (VFNS) @, 0|- The VFNS is an approach alternative to the traditional fixed flavor 
number scheme (FFNS) where only light degrees of freedom (it, d, s and g) are considered as active. It is well known 
that a heavy quark production cross section contains potentially large logarithms of the type a s In (Q 2 /to 2 ) whose 
contribution dominates at high energies, Q 2 — > oo. Within the VFNS, these mass logarithms are resummed through 
the all orders into a heavy quark density which evolves with Q 2 according to the standard DGLAP [1, [ij [l(| evolution 
equation. Hence the VFN schemes introduce the parton distribution functions (PDFs) for the heavy quarks and 
change the number of active flavors by one unit when a heavy quark threshold is crossed. We can say that the charm 
density arises within the VFNS perturbatively via the g — > cc evolution and will call it the perturbative IC. 

Presently, both perturbative and nonperturbative IC are widely used for a phcnomcnological description of available 
data. (A recent review of the theory and experimental constraints on the charm quark distribution can be found in 
Refs. [ljlllll]- See also AppendixOin the present paper). In particular, practically all the recent versions of the CTEQ 
[l3[ and MRST [l4[ sets of PDFs are based on the VFN schemes and contain a charm density. At the same time, the 
key question remains open: How to measure the intrinsic charm content of the proton? The basic theoretical problem 
is that radiative corrections to the fixed order predictions for the production cross sections are large. In particular, 
the next-to- leading order (NLO) corrections increase the leading order (LO) results for both charm and bottom 
production cross sections by approximately a factor of two at energies of the fixed target experiments. Moreover, soft 
gluon rcsummation of the threshold Sudakov logarithms indicates that higher-order contributions are also essential. 
(For a review see Refs. [H|,[l6j]). On the other hand, perturbative instability leads to a high sensitivity of the theoretical 
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calculations to standard uncertainties in the input QCD parameters. For this reason, it is difficult to compare pQCD 
results for spin-averaged cross sections with experimental data directly, without additional assumptions. The total 
uncertainties associated with the unknown values of the heavy quark mass, m, the factorization and renormalization 
scales, fip and fiR, Aqcd and the PDFs are so large that one can only estimate the order of magnitude of the pQCD 
predictions for production cross sections [lT], Ell • 

Since production cross sections are not perturbatively stable, it is of special interest to study those observables that 
are well-defined in pQCD. A nontrivial example of such an observable was proposed in Refs. [HI, [2(| HH [HJ where 
the azimuthal cos2ip asymmetry in heavy quark photo- and leptoproduction has been analyzed 1 . In particular, the 
Born level results have been considered [19l| and the NLO soft-gluon corrections to the basic mechanism, photon- 
gluon fusion (GF), have been calculated [lO, Hl|- It was shown that, contrary to the production cross sections, 
the azimuthal asymmetry in heavy flavor photo- and leptoproduction is quantitatively well defined in pQCD: the 
contribution of the dominant GF mechanism to the asymmetry is stable, both parametrically and perturbatively. 
Therefore, measurements of this asymmetry would provide an ideal test of pQCD. As was shown in Ref. [22T |. the 
azimuthal asymmetry in open charm photoproduction could have been measured with an accuracy of about ten percent 
in the approved E160/E161 experiments at SLAC [23| using the inclusive spectra of secondary (decay) leptons. 

In the present paper we study the IC contribution to the azimuthal asymmetries in heavy quark leptoproduction: 

l(£) + N(p) -> 1(1 -q) + Q[p Q ) + X[Q]( PX ). (1) 

Neglecting the contribution of Z— boson as well as the target mass effects, the cross section of the reaction (fl]) for 
unpolarizcd initial states may be written as 

d 3 a m a em 1 y 2 



a T (x,Q 2 ) + ea L (x,Q 2 ) + ea A (x, Q 2 ) cos 2ip + 2- v /e(l + e)ai(x, Q 2 ) cos tp , (2) 



dxdQ 2 d^ (2tt) 2 xQ 2 1 - e 
The quantity e measures the degree of the longitudinal polarization of the virtual photon in the Breit frame [24 

^ l + (l-2/) 2 ' 



(3) 



and the kinematic variables are defined by 

S = (£ + p) 2 , Q 2 = -g 2 , 



Q 2 



2p-q' 



P'Q rvi a 4m 2 

V= — j, Q 2 = xyS, P=- 3 ~. 4 
p ■ I S 

The cross sections (Xj (i — T, L, A, I) in Eq. ^ are related to the structure functions Fi(x, Q 2 ) as follows: 



F t (x,Q 2 ) = f cn(x,Q 2 ), (i = T,L,A,I) 



F 2 (x,Q 2 ) - — ^«7 2 (x,Q 2 ), (5) 



4tt 2 



where F 2 — 2x(Ft+Fl) and a 2 = ctt + vl- In Eq. @, ctt (ox) is the usual 7*^V cross section describing heavy quark 
production by a transverse (longitudinal) virtual photon. The third cross section, a a, comes about from interference 
between transverse states and is resp onsible for the cos 2tp asymmetry which occurs in real photoproduction using 
linearly polarized photons [l9ll20l[22| . The fourth cross section, crj, originates from interference between longitudinal 
and transverse components |24l| In the nucleon rest frame, the azimuth ip is the angle between the lepton scattering 
plane and the heavy quark production plane, defined by the exchanged photon and the detected quark Q (see Fig.[TJ. 
The covariant definition of tp is 



r-n . Q 2 ^l/x 2 +4m 2 N /Q 2 

cos <p = —j=^-j= , sin tp — — = n ■ t, (6) 

\j — r 2 \J — n 2 2\/—r 2 \/—n 2 

r» = e^ a Pp„q a e Pl n» = e^q uPa p Q0 . (7) 



1 The well-known examples are the shapes of differential cross sections of heavy flavor production which are sufficiently stable under 
radiative corrections. 




FIG. 1: Definition of the azimuthal angle tp in the nucleon rest frame. 



In Eqs. f4| and ((6j> , m and are the masses of the heavy quark and the target, respectively. Usually, the azimuthal 
asymmetry associated with the cos2tp distribution, A2 V (p,x,Q 2 ), is defined by 

\2\ o/ o u _ d3<J iN(<P = 0) + d 3 a lN (tp = tt) - 2d 3 ai N (ip = vr/2) 



A 2v (p,a;,(5 2 ) = 2(cos2^)(p,2:,Q 2 ) = 



d 3 ai N (ip = 0) + d 3 er ;A r(^ = 7r) + 2d 3 cr;Ar(v3 = n/2) 



eo A (x,Q 2 ) _ A{xQ2) e + eR{x,C?) 



a T (x, Q 2 ) + ea L {x,Q 2 ) ' l + eR(x,Q 

(n. 

dxdQ 2 dtp 



where d 3 aiN{<p) = ^ — (p, x, Q 2 , tp) and the mean value of cosntp is 



2 - d 3 



/ dtp cos ntp—^-^—(p,x,Q 2 ,tp) 



(cos ntp)(p 7 x,Q 2 ) = dxd Q d 1 _ 

[A d ° lN I 

I (f dxdQ 2 dtp { 

In Eq. ((U, the quantities R(x,Q 2 ) and A(x,Q 2 ) are defined as 



R(x,Q 2 ) = —(x,Q 2 ) = ^(x,Q 2 ), (10) 
A(x,Q 2 ) = ^A( x ,Q 2 ) = 2x^(x,Q 2 ). (11) 

<J2 t<2 

Likewise, we can define the azimuthal asymmetry associated with the cos tp distribution, A v (p, x, Q 2 )'. 

2d 3 a lN (tp = 0) - 2d 3 cr w (^ = tt) 



A v (p,x,Q 2 ) = 2(cos tp}(p, x, Q 2 ) = 



d 3 (Ji N (tp = 0) + d 3 <r lN (tp = tt) + 2d 3 ai N (tp = w/2) 



a T (x,Q 2 ) + ea L {x,Q 2 ) l + eR(x,Q 2 ) 

where 

A x {x, Q 2 ) = 2V2 — (x, Q 2 ) = AV2x^-(x, Q 2 ). (13) 

C2 ^2 

Remember that y <C 1 in most of the experimentally reachable kinematic range. Taking also into account that 
e= l + 0(y 2 ), we find: 

A 2v ( P ,x,Q 2 ) = A(x,Q 2 )+0(y 2 ), A v (p,x,Q 2 ) = Aj{x,Q 2 ) + 0(y 2 ). (14) 
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So, like the <J2(x,Q 2 ) cross section in the (^-independent case, it is the parameters A(x,Q 2 ) and Aj(x, Q 2 ) that can 
effectively be measured in the azimuth-dependent production. 

In this paper we concentrate on the azimuthal asymmetry A(x,Q 2 ) associated with the cos 2<p-distribution. We 
have calculated the IC contribution to the asymmetry which is described at the parton level by the photon-quark 
scattering (QS) mechanism given in Fig. [H Our main result can be formulated as follows: 

• Contrary to the basic GF component, the IC mechanism is practically cos 2</?-independent. This is due to the 
fact that the QS contribution to the cta(x, Q 2 ) cross section is absent (for the kinematic reason) at LO and is 
negligibly small (of the order of 1%) at NLO. 

As to the ^-independent cross sections, our parton level calculations have been compared with the previous results for 
the IC contribution to (Jzix, Q 2 ) and cjl{x, Q 2 ) presented in Refs. [25l[26|. Apart from two trivial misprints uncovered 
in [25| for <tl{x, Q 2 ), a complete agreement between all the considered results is found. 

Since the GF and QS mechanisms have strongly different cos 2^-distributions, we investigate the possibility to 
discriminate between their contributions using the azimuthal asymmetry A(x,Q 2 ). We analyze separately the non- 
perturbative IC in the framework of the FFNS and the perturbative IC within the VFNS. 

The following properties of the nonperturbative IC contribution to the azimuthal asymmetry within the FFNS are 
found: 

• The nonperturbative IC is practically invisible at low x, but affects essentially the GF predictions at large x. 
The dominance of the cos 2^-independent IC component at large x leads to a more rapid (in comparison with 
the GF predictions) decreasing of A(x, Q 2 ) with growth of x. 

• Contrary to the production cross sections, the cos 2<p asymmetry in charm azimuthal distributions is practically 
insensitive to radiative corrections at Q 2 ~ m 2 . Perturbative stability of the combined GF+QS result for 
A(x, Q 2 ) is mainly due to the cancellation of large NLO corrections in Eq. (fTTj) . 

• pQCD predictions for the cos2c/? asymmetry are parametrically stable; the GF+QS contribution to A(x,Q 2 ) is 
practically insensitive to most of the standard uncertainties in the QCD input parameters: (j,r, /i^, Aqcd and 
PDFs. 

• Nonperturbative corrections to the charm azimuthal asymmetry due to the gluon transverse motion in the target 
are of the order of 20% at Q 2 < m 2 and rapidly vanish at Q 2 > m 2 . 

We conclude that the contributions of both GF and IC components to the cos 2ip asymmetry in charm leptoproduction 
are quantitatively well defined in the FFNS: they are stable, both parametrically and perturbatively, and insensitive 
(at Q 2 > m 2 ) to the gluon transverse motion in the proton. At large Bjorken x, the A(x, Q 2 ) asymmetry could be a 
sensitive probe of the nonperturbative IC. 

The perturbative IC has been considered within the VFNS proposed in Refs. [f| Q. The following features of the 
azimuthal asymmetry should be emphasized: 

• Contrary to the nonperturbative IC component, the perturbative one is significant practically at all values of 
Bjorken x and Q 2 > m 2 . 

• The charm densities of the recent CTEQ and MRST sets of PDFs lead to a sizeable reduction (by about 1/3) 
of the GF predictions for the cos 2ip asymmetry. 

We conclude that impact of the perturbative IC on the cos2</? asymmetry is sizeable in the whole region of x and, for 
this reason, can easily be detected. 

Concerning the experimental aspects, azimuthal asymmetries in charm leptoproduction can, in principle, be mea- 
sured in the COMPASS experiment at CERN, as well as in future studies at the proposed eRHIC [23, [28[ and LHeC 
[211] colliders at BNL and CERN, correspondingly. 

The paper is organized as follows. In Section |TT] we analyze the QS and GF parton level predictions for the in- 
dependent charm leptoproduction in the single-particle inclusive kinematics. In particular, we discuss our results for 
the NLO QS cross sections and compare them with available calculations. Hadron level predictions for A(x, Q 2 ) are 
given in Section [mi We consider the IC contributions to the asymmetry within the FFNS and VFNS in a wide region 
of x and Q 2 . Some details of our calculations of the QS cross sections are presented in Appendix An overview of 
the soft-gluon resummation for the photon-gluon fusion mechanism is given in Appendix [Bl Some experimental facts 
in favor of the nonperturbative IC are briefly listed in Appendix \C\ 
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FIG. 2: The LO (a) and NLO (b and c) photon-quark scattering diagrams. 



II. PARTONIC CROSS SECTIONS 
A. Quark Scattering Mechanism 

The momentum assignment of the deep inelastic lepton-quark scattering will be denoted as 

m + Q(k Q ) -> l(£ -q) + Q( PQ ) + X(p x ). (15) 
Taking into account the target mass effects, the corresponding partonic cross section can be written as follows [3C| 



1 oiq ot em y 2 \/T+ 4Az 2 



dzdQMip (2tt) 2 zQ 2 1 - e 



& 2 ,q{z,X) - (1 - e)a L . Q (z, A) + ia AtQ {z, A) cos 2^ + 2yje(l +i)& I , Q (z, A) cos^ . 

(16) 



In Eq. (|16p . we use the following definition of partonic kinematic variables: 

Q ■ kg Q 2 m 2 _ 

£■ k Q 2q-k Q Q 2 

In the massive case, the (virtual) photon polarization parameter, £, has the form [3 

2(1 -y- Xz 2 y 2 ) 



(18) 



1 + (1 - y) 2 + 2Xz 2 y 2 ' 
At leading order, 0(a em ), the only quark scattering subprocess is 

l*(q) + Q(k Q )^Q(p Q ). (19) 
The 7*Q cross sections, a^L (A; = 2, L, A, I), corresponding to the Born diagram (see Fig. [2^,) are: 



4°q(z,X) = v B {z)\Tl + 4Az 2 5(1 -z), 

ai 0) Q (z,A) = ^(z) 4AZ ' ,5(1 -z), (20) 

V 1 + 4:\Z Z 

a%(z,X) = &fUz,X) = 0, 
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with 



<tb(z) 



(2n) 2 e 2 Q a em 

Q 2 



(21) 



where eQ is the quark charge in units of electromagnetic coupling constant. 

To take into account the NLO 0(a em a s ) contributions, one needs to calculate the virtual corrections to the Born 
process (given in Fig. [2k) as well as the real gluon emission (see Fig.[2h>): 



i*(q) + Q(k Q ) ^ Q(pq) + g{p g 



(22) 



The NLO (^-dependent cross sections, cr^Q an< ^ ^iQi are described by the real gluon emission only. Corresponding 

contributions are free of any type of singularities and the quantities and o^ q can be calculated directly in four 
dimensions. 

In the ^-independent case, a^Q and &lq, we also work in four dimensions. The virtual contribution (Fig. 
contains ultraviolet (UV) singularity that is removed using the on-mass-shell regularization scheme. In particular, we 
calculate the absorptive part of the Feynman diagram which has no UV divergences. The real part is then obtained 
by using the appropriate dispersion relations. As to the infrared (IR) singularity, it is regularized with the help of 
an infinitesimal gluon mass. This IR divergence is cancelled when we add the bremsstrahlung contribution (Fig. [2Jd). 
Some details of our calculations are given in Appendix [X] 

The final (real+ virtual) results for 7*Q cross sections can be cast into the following form: 

a£l(z, A) = ^C F a B (l)VT+4XS(l - z){-2 + 4 In A - VTTlA lnr + - + 2A [2Li 2 (r 2 ) + 4Li 2 (-r) 



VI + 4A 

+31n 2 (r) - 41nr + 41nrln(l + 4A) - 21nrlnA j 



-^ CfMz) (1+4\z>)V> {[1-(1- AM 



1 - 3z - 4z 2 + 6z 3 + 8z 4 - 8z D 
+6Az (3 - 18z + 13z 2 + 10z 3 - 8z 4 ) 

+4A 2 z 2 (8 - 77z + 65z 2 - 2z 3 ) 
+16A 3 z 3 (l - 21z + 12z 2 ) - 128A 4 z 5 



2 lnD(z, A) 
' VI + 4Az 2 



-(1 + Z + 2Z 2 



2z A 



2\z (2 - \\z - llz 2 ) + 8A 2 z 2 (1 - 9z) 



8(1 + 4A) 2 z 4 4(1 + 2A)(1 + 4A) 2 z 4 bxD(z, A) 



(!-*)_ 



Vl + 4Az 2 



(1-*)h 



(23) 



^ Q (z,X) = -C F a B (l) 7 == 



5(1 -z){ 



-2 + 41nA - 



4A 



In / 



1 + 2A 



2Li 2 (r 2 )+4Li 2 (-r) 



, as n - I \ 1 J z 

+ -C fCTb (z )(i + 4Az2)3/2 | [i _ (i _ a ^ 



vT + 4A VI + 4A 

+31n 2 (r) -41nr + 41nrln(l + 4A) - 21nrlnA } 

(1 - z) 2 - Az (13 - 19z - 2z 2 + 8z 3 ) 



-2A 2 z 2 (31-39z + 8z 2 ) 

-8A 3 z 3 (10 - 7z) - 32A 4 z 4 

2Az 2 lnL>(z,A) r , , 

v ' ; 3 + 3z + 16Az 

Vl + 4Az 2 L J 

8A(1 + 4A)z 4 4A(1 + 2A)(1 + 4A)z 4 In D(z, A) 



(l-z)_ 



VI +4Az 2 



(l-z)_ 



(24) 



°A,qM = ^ C ^^(i + 4Az 2 )3/4[l-(l-A)z] 



[l + 2A(4-3z) + 8A 2 z] + 



2AlnD(z, A) 
VI + 4Az 2 



[2 + z + 4Xz]\, (25) 
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8V2 



C F a B (z) 



(l + 4Az 2 ) 2 [i-(i-X)z 



3/2 



-(1 - z)(l + 2z) - 4Az (10 - Wz - z 2 + 2z 3 ) 



-8\ 2 z 2 (25 - 29z + 8z 2 ) - 96A 3 z 3 (3 - 2z) - 128AV (26) 
+8 yj\z [1 - (1 - Xjz] [1 - z 2 + Az(13 - llz) + 4AV(7 - Az) + 16A 3 z 3 ] 



In Eqs. (|23H26p , C F = (N 2 - 1)/(2N C ), where iV c is number of colors, while 



D(z,X) = 



1 + 2Az - VI + 4Az 2 



1 + 2Az + VI + 4Az 2 
The so-called "plus" distributions are defined by 



r = vAD(T~T7A) = 



V1 + 4A- 1 



V1 + 4A+ 1 



(27) 



1 

[g(z)] + = g(z) - 5(1 - z) J dCs(C)- 



(28) 



For any sufficiently regular test function h(z), Eq. I|28p gives 

l r , -i l 

dz h(z) 



ln fe (l- z) 



1- z 



, ln fe (l-z) r , , . , In +1 (1 — a) 

dz — r \h(z) - h(l)} + h(l)- 



1 - z 



k + l 



(29) 



To perform a numerical investigation of the inclusive partonic cross sections, 07c, q (k — T, L, A, I), it is convenient 



to introduce the dimensionless coefficient functions 



>.0 

-fc,Q > 



2 /2\oo n / 2 



(30) 



where fi is a factorization scale (we use (J, = fiF = M_r) and the variable 77 measures the distance to the partonic 
threshold: 



s 1 — z 

^ = ~~ 2 ~ 1 = ~ 5 — ' 

m z Az 



S = (q + kg) 2 



(31) 



Our analysis of the quantity c^'q (t}, A) is given in Fig. [31 One can see that c^'q is negative at low Q 2 (A -1 < 1) 

and positive at high Q 2 (A -1 > 20). For the intermediate values of Q 2 , c^'q(ji, A) is an alternating function of 77. 

Our results for the coefficient function Cj °q (77, A) at several values of A are presented in Fig. [3] It is seen that Cjq 
is negative at all values of 77 and A. Note also the threshold behavior of the coefficient function: 

4° Q 0) (77 - 0, A) = -V2n 2 C FJ ^ + 0{ V ). (32) 

This quantity takes its minimum value at A m = 1/4: Cj q (77 = 0, A m ) = —it 2 Gf/ (2y/2). 



B. Comparison with Available Results 

For the first time, the NLO 0(a em a s ) corrections to the ^-independent IC contribution have been calculated a 
long time ago by Hoffmann and Moore (HM) t 25]. However, authors of Ref. [25| don't give explicitly their definition 
of the partonic cross sections that leads to a confusion in interpretation of the original HM results. To clarify the 
situation, we need first to derive the relation between the lepton-quark DIS cross section, &<jiq, and the partonic cross 
sections, and a^ L \ used in 25]. Using Eqs. (C.l) and (C.5) in Ref. [25| . one can express the HM tensor cr^ in 
terms of "our" cross sections o"2,q and ctl,q defined by Eq. (|16p in the present paper. Comparing the obtained results 
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FIG. 3: c^'q (r), A) and Cj°q' (77, A) coefficient functions at several values of A. 



with the corresponding definition of via the HM cross sections and (given by Eqs. (C.16) and (C.17) in 
Ref. [Hj]), we find that 



02,qO,A) 
°l,q(z, A) 



0\b 



(z)Vl + 4AzV 2 >(z,A), 



25/3 (z) 
VI + 4Az 2 



/ L )(z,A) + 2AzV 2 )(^A) 



(33) 
(34) 



Now we are able to compare our results with original HM ones. It is easy to see that the LO cross sections (defined 
by Eqs. (37) in and Eqs. (f20|) in our paper) obey both above identities. Comparing with each other the quantities 

crp' ) and a^Q (given by Eq. (51) in [IE] and Eq. (f23|) in this paper, respectively), we find that identity (l33|) is satisfied 
at NLO too. The situation with longitudinal cross sections is more complicated. We have uncovered two misprints in 
the NLO expression for given by Eq. (52) in j25j. First, the r.h.s. of this Eq. must be multiplied by z. Second, 
the sign in front of the last term (proportional to 6(1 — z)) in Eq. (52) in Ref. 25] must be changed 2 . Taking into 
account these typos, we find that relation ([34| holds at NLO as well. So, our calculations of &2,Q and <Jl,q agree with 
the HM results. 

Recently, the heavy quark initiated contributions to the (/^-independent DIS structure functions, F2 and Fl, have 
been calculated by Kretzer and Schienbein (KS) [26J. The final KS results are expressed in terms of the parton level 
structure functions Hf and iff- Using the definition of Hf and iff given by Eqs. (7, 8) in Ref. [26(, we obtain that 



&t,q(z, A) 



a s d B {z) HfjCX) 
2tt VI + 4A Vl + 4Az 2 ' 



&2,q(z, A) = -^-vb(z) 
Zir 



1 + 4A 
1 + 4Az 2 



A), (35) 



where ctt,q = <J2,Q — &l,q and <tl,q are defined by Eq. (flUj) in our paper and £' = z (l + VI + 4A) /(l + VT + 4Az 2 ) . 
To test identities (I33t . one needs only to rewrite the NLO expressions for the functions Hf(£', A) and H^iC , A) (given 
in Appendix C in Ref. [26]) in terms of variables z and A. Our analysis shows that relations (|35[) hold at both LO 
and NLO. Hence we coincide with the KS predictions for the 7*Q cross sections. 

However, we disagree with the conclusion of Ref. [26| that there are errors in the NLO expression for given in 
Ref. [25[ 3 . As explained above, a correct interpretation of the quantities and used in [25| leads to a complete 
agreement between the HM, KS and our results for ^-independent cross sections. 

As to the (^-dependent DIS, pQCD predictions for the 7*Q cross sections aA,Q{z,X) and 07,(3(2:, A) in the case 
of arbitrary values of m and Q 2 are not, to our knowledge, available in the literature. For this reason, we have 
performed several cross checks of our results against well known calculations in two limits: m 2 — > and Q 2 — * 0. 
In particular, in the chiral limit, we reproduce the original results of Georgi and Politzer (32j and Mendez [33[ for 



2 Note that this term originates from virtual corrections and the virtual part of the longitudinal cross section given by Eq. (39) in Ref. [25I ] 
also has wrong sign. See Appendix^for more details. 

3 In detail, the KS point of view on the HM results is presented in PhD thesis [311 ]. pp. 158-160. 
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fJ 



FIG. 4: The LO photon-gluon fusion diagrams. 

&i t Q(z,X — > 0) and &a,q(z, A — > 0). In the case of Q 2 — > 0, our predictions for <72,q(s,Q 2 — * 0) and &a,q{s,Q^-^ 0) 
given by Eqs. (|23I25[) reduce to the QED textbook results for the Compton scattering of polarized photons [34| . 

C. Photon-Gluon Fusion 

The gluon fusion component of the semi-inclusive DIS is the following parton level interaction: 

l(£) + g(k g ) -> l(£ -q) + Q(p Q ) + X[Q](p x ). (36) 

Corresponding lepton-gluon cross section, d<r/ g , has the following decomposition in terms of the helicity j*g cross 
sections: 

3~, ™ 1 -,2 



d J <Ti g a em 1 y 



dzdQ 2 dtp (2tt) 2 zQ 2 1 - e 



&2, g (z,\) - (1 - e)a L ,g(z, A) + eaA,g(z, A) cos2<£> + 2a/£(1 + e)ai, g (z, A) cosy , (37) 



where the quantity e is defined by Eq. |3]) with y — (q ■ kg) j (£ ■ k g ). 

At LO, 0(a em a s ), the only gluon fusion subprocess responsible for heavy flavor production is 

r(q)+g(k g )^Q(PQ) + Q( P Q). (38) 

The 7*g cross sections, djp (k — 2, L. A, I), corresponding to the Born diagrams given in Fig.0]have the form [35l. l36l] : 

A ) = ^b{z){ [(1 - zf + z 2 + 4Az(l - 3z) - 8AV] In - [1 + 4z(l - z)(A - 2)] 

4% A) = ^ B (^{-2Azlni±|i + (l-z)&}, 

4°i(^A) = ^a B (z)z{2A[l-2z(l + A)]lni±|i + (l-2A)(l-z)/3 z }, (39) 

*S(*>*) = °. 

where <Jb(z) is defined by Eq. (121|) and the following notations are used: 



x m2 a L 4AT 

2= 2^- x= q^ ^=r-— (40) 

Note that the cosy dependence vanishes in the GF mechanism due to the Q <-> Q symmetry which, at leading order, 
requires invariance under <p — > <p + it [37j ■ 

As to the NLO results, presently, only <p- independent quantities a^ g and er^ 1 ^ are known exactly (38|. For this 
reason, we will use in our analysis the so-called soft-gluon approximation for the NLO j*g cross sections (see Appendix 
|B|) . As shown in Refs. [2(| HJ, H^] , at energies not so far from the production threshold, the soft-gluon radiation is 
the dominant perturbative mechanism in the j*g interactions. 
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III. HADRON LEVEL RESULTS 



A. Fixed Flavor Number Scheme and Nonperturbative Intrinsic Charm 

In the fixed flavor number scheme 4 , the wave function of the proton consists of light quarks u,d,s and gluons 
g. Heavy flavor production in DIS is dominated by the gluon fusion mechanism. Corresponding hadron level cross 
sections, dk,GF(x,X), have the form 

1 

a k7GF (x,X) = Jdzg(z,n F )a k , g {x/z,\,fi F ), (k = 2,L,A,I), (41) 

x 

X = x(l + 4X), (42) 

where g(z,/j,p) describes gluon density in the proton evaluated at a factorization scale //p. The lowest order GF 
cross sections, (k = 2,L,A,I), are given by Eqs. (|39[) . The NLO results, d^ g , to the next-to-leading logarithmic 
accuracy are presented in Appendix [B] 

We neglect the j*q(q) fusion subprocesses. This is justified as their contributions to heavy quark leptoproduction 
vanish at LO and are small at NLO [38]. 

In the FFNS, the intrinsic heavy flavor component of the proton wave function is generated by gg — > QQ fluctuations 
where the gluons are coupled to different valence quarks. In the present paper, this component is referred to as the 
nonperturbative intrinsic charm (bottom). The probability of the corresponding five-quark Fock state, \uudQQ}, is 
of higher twist since it scales as Aq CD /m 2 @. However, since all of the quarks tend to travel coherently at same 
rapidity in the \uudQQ^ bound state, the heaviest constituents carry the largest momentum fraction. For this reason, 
the heavy flavor distribution function has a more "hard" z-behavior than the light parton densities. Since all of the 
densities vanish at z — ► 1, the hardest PDF becomes dominant at sufficiently large z independently of normalization. 

Convolution of PDFs with partonic cross sections does not violate this observation. In particular, assuming a gluon 
density g(z) ~ (1 — z) n (where n = 3 — 5), we obtain that the LO GF contribution to F2 scales as (1 — ^-) n + 3 / 2 
at x ~ * 1) where x iS defined by Eq. ((4"2"|) . In the case of Hoffman and Moore charm density (see below), the LO 
IC contribution is proportional to (1 — %) at % — » 1. It is easy to see that, independently of normalizations, the IC 
contribution to be dominate over the more "soft" GF component at large enough x. 

For the first time, the intrinsic charm momentum distribution in the five-quark state \uudcc) was derived by Brodsky, 
Hoyer, Peterson and Sakai (BHPS) in the framework of a light-cone model [3, Neglecting the transverse motion of 
constituents, they have obtained in the heavy quark limit that 

c(z) = ^-z 2 \6z(l + z)ln z + (1 - z)(l + lOz + z 2 )] , (43) 
6 

where N$ = 36 corresponds to a 1% probability for IC in the nucleon: J Q c(z)dz = 0.01. 

Hoffmann and Moore (HM) [25[ incorporated mass effects in the BHPS approach. They first introduced a mass 
scaling variable £, 



2ax 1 + V1 + 4A , AA . 
f = ■ a= , (44) 

where A at = m 2 N j Q 2 . To provide correct threshold behavior of the charm density, the constraint £ < 7 < 1 was 
imposed where 

2ai 1 (45) 



1 + V1 + 4Ajvx 2 1 + 4A-Aat 

Resulting charm distribution function, c(£, 7), has the following form in the HM approach 



Mi.-»=r (0 ~7 C(7) ' (46) 
0, £> 7 



4 This approach is sometimes referred to as the fixed-order perturbation theory (FOPT). 
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with c(£) defined by Eq. (|43p . Corresponding hadron level cross sections for the (c + c) production, (Jk,Qs{x 1 A), due 
to the heavy quark scattering (QS) mechanism, are 



<?k,Qs{x, A) = 



d2 



;c+0*>7) ^fc,c(£/z, A) , 



(fc = 2,L,A,7), 



(47) 



where the charm density c+(z,7) = 0(2:, 7) + 0(2,7). The LO and NLO expressions for the partonic cross sections 
<Jk, c (z 7 A) are given by Eqs. (|2TJ|) and (|23H26j) . respectively. 

Note also that, in the FFNS, the full cross section for the charm production, <Tk(x, A), is simply a sum of the GF 
and IC components: 



a k (x, A) = a k , G F{x, A) + a k:QS (x, A), 



(k = 2,L,A,I). 



(48) 



Let us discuss the FFNS predictions for the hadron level asymmetry parameter A(x, Q 2 ) defined by Eq. (jTTJ) . First 
we consider the ratio (cr^QS / (t 2 > qs > ) ( x > i- e -' tne mere IC contribution to A(x,Q 2 ). In Fig. El the x-behavior of 
this quantity at various values of A is presented. One can see that the ratio (ct^qs /°2°qs^) ( x > i s negligibly small 
(of the order of 1%) practically at all values of Q 2 > m 2 . Note that this fact is independent of the charm density 
we use (see, for instance, the right panel in Fig. O, but is only due to the smallness of the partonic cross section 

a^ c (z, A) [30( | . So, the quantity ct j 4,qs(^,A) is exactly zero at LO 5 and negligibly small at NLO. This implies that 
the IC contribution has no practically cos 293-dependence and, for this reason, we will neglect both &a,c(z, A) and 
&a,qs(x, A) cross sections in our further analysis. 

Fig.E]shows A(x, A) as a function of x for several values of variable A: A -1 = 1, 4, 10, 20, 50 and 100. We display both 
LO and NLO predictions of the GF mechanism as well as the analogous results of the combined GF+QS contribution. 
The azimuthal asymmetry due to the mere LO GF component is given by solid line. The NLO GF predictions are 
plotted by dashed line. The LO and NLO results of the total GF+QS contribution are given by dash-dotted and long- 
dashed lines, respectively. In our calculations, the CTEQ5M [4(| parametrization of the gluon distribution function is 
used and a 1% probability for IC in the nucleon is assumed. Throughout this paper, the value fiF — I^r — y fn 2 + Q 2 
for both factorization and renormalization scales is chosen. In accordance with the CTEQ5M parametrization, we 
use m c = 1.3 GeV and A 4 = 326 MeV H}. 

One can see from Fig. [S] the following basic features of the azimuthal asymmetry, A(x, A), within the FFNS. First, 
as expected, the nonperturbative IC contribution is practically invisible at low x, but affects essentially the GF 
predictions at large x. Since, contrary to the GF mechanism, the QS component is practically cos 2tp- independent, 



5 Remember that the LO quantity a ? (z, A) vanishes for the kinematic reason, see Eqs. (1206 . 
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the dominance of the IC contribution at large x leads to a more rapid (in comparison with the GF predictions) 
decreasing of A(x, A) with growth of x. 

The most remarkable property of the azimuthal asymmetry is its perturbative stability. In Refs. [I^jHHj the NLO 
soft-gluon corrections to the GF predictions for the cos 2tp asymmetry in heavy quark photo- and leptoproduction was 
calculated. It was shown that, contrary to the production cross sections, the quantity A(x, A) is practically insensitive 
to soft radiation. One can see from Fig. [6] in the present paper that the NLO corrections to the LO GF predictions 
for A{x, A) are about few percent at not large x. This implies that large soft-gluon corrections to <J^gf an< ^ ^gf 
(increasing both cross sections by a factor of two) cancel each other in the ratio {pagf/ °2GF) fa' •V) wu -h a good 
accuracy. In terms of so-called K -factors, Kk(x, A) = (^jA^'A (t^°' > )(x, A) for k = 2, L, A, I, perturbative stability of 
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the GF predictions for A(x, A) is provided by the fact that the corresponding AT-factors are approximately the same 
at not large x: Ka,gf{%, A) w K2,gf(x, A). 

Comparing with each other the dash-dotted and long-dashed curves in Fig. [51 we see that the NLO corrections to the 
combined GF+QS result for A(x, A) are also small. In this case, three reasons are responsible for the closeness of the 
LO and NLO predictions. At small x, where the nonperturbative IC contribution is negligible, perturbative stability 
of the asymmetry is provided by the GF component. In the large- x region, where the IC mechanism dominates, the 
azimuthal asymmetry rapidly vanishes with growth of x at both LO and NLO because the QS component is practically 
cos2<p-independent, <?a\{x, A) « cr^t c (#j A) = 6 . At intermediate values of x, where both mechanisms are essential, 
perturbative stability of A{x, A) is due to the similarity of the corresponding if -factors: K2,gf{%, A) ~ K2,qs{%, A) 7 . 

Another remarkable property of the azimuthal asymmetry closely related to fast perturbative convergence is its 
parametric stability 8 . The analysis of Refs. [l^, HH shows that the GF predictions for the cos2tp asymmetry are less 
sensitive to standard uncertainties in the QCD input parameters (to, [ir, fip, Aqcd and PDFs) than the corresponding 
ones for the production cross sections. We have verified that the same situation takes also place for the combined 
GF+QS results. 

Let us discuss how the GF predictions for the azimuthal asymmetry are affected by nonperturbative contributions 
due to the intrinsic transverse motion of the gluon in the target. Because of the relatively low c-quark mass, these 
contributions are especially important in the description of the cross sections for charmed particle production. 

To introduce fey degrees of freedom, k g ~ £p + fey, one extends the integral over the parton distribution function 
in Eq. fH]) to fey-space, 

d(g(C,VF)^d(d 2 k T f(k T )g(<; if i F ). (49) 
The transverse momentum distribution, /(fcr), is usually taken to be a Gaussian: 

f w = ~^mr- (50) 

In practice, an analytic treatment of fey effects is usually used. According to Ref. [ll|, the fey-smeared differential 
cross section of the process ([T]) is a two-dimensional convolution: 

dxdQMpQTdtp = / d2kT \(k* T ) dxdQ*dp Q Tdip(? QT - ¥ T ) • (51) 

The factor ~ in front of fey in the r.h.s. of Eq. (fSTj) reflects the fact that the heavy quark carries away about one half 
of the initial energy in the reaction |T]). 

Values of the fey-kick corrections to the LO GF predictions for the cos 2<p asymmetry in the charm production 
are shown in Fig. [6] by dotted curves. Calculating the fey-kick effect we use (fe T ) = 0.5 GeV 2 . One can see that 
fey-smearing for A(x,Q 2 ) is about 20-25% in the region of low Q 2 < m 2 and rapidly decreases at high Q 2 . 

In Fig. [7J the dependence of the asymmetry A(x, A) on the nonperturbative intrinsic charm content of the proton is 
presented. We plot the LO predictions for A(x. A) as a function of x for several values of the variable A and quantity 
P c = j Q c(z)dz describing a probability for IC in the nucleon. Dash-dotted curves describe the GF( L °)+QS( LO) 
contributions with P c = 5%, 1%, 0.1% and 0.01%. Solid lines correspond to the case when P c — 0. Comparing with 
each other Figs. [6jand[7l one can see that even a 0.1% contribution of the nonperturbative IC to the proton wave 
function could be extracted from the cos 2ip asymmetry at large enough Bjorken x. 



B. Variable Flavor Number Scheme and Perturbative Intrinsic Charm 

One can see from Eqs. (|3"9")l that the GF cross section a^^z^X) contains potentially large logarithm, ln(Q 2 /m 2 ). 
The same situation takes also place for the QS cross section ^^(z, A) given by Eq. (f23|) . At high energies, Q 2 — > 



6 Although the ratio {A^ NL °^ /A^ L °^){x,\) is sizeable at sufficiently large x, the absolute values of the quantities A^ L °^(x,X) and 
A( nlo ^(x, A) become so small that it seems reasonable to consider the asymmetry as equally negligible at both LO and NLO and treat 
the predictions as perturbatively stable. 

7 Note however that this similarity takes only place at intermediate values of x where both GF and QS components are essential. In the 
low- and large- a; regions, the factors -R'2,Gf( 3; ! -*0 an d K 2,Qs{ x > A) are strongly different. 

8 Of course, parametric stability of the fixed order results does not imply a fast convergence of the corresponding series. However, a fast 
convergent series must be parametrically stable. In particular, it must be fin- and /ip -independent. 
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oo, the terms of the form a s ln(Q 2 /m 2 ) dominate the production cross sections. To improve the convergence of 
the perturbative series at high energies, the so-called variable flavor number schemes (VFNS) have been proposed. 
Originally, this approach was formulated by Aivazis, Collins, Olness and Tung (ACOT) [f| |42j • 

In the VFNS, the mass logarithms of the type a™ ln n (Q 2 /m 2 ) are resummed via the renormalization group equations. 
In practice, the resummation procedure consists of two steps. First, the mass logarithms have to be subtracted from 
the fixed order predictions for the partonic cross sections in such a way that in the asymptotic limit Q 2 — > oo the well 
known massless MS coefficient functions are recovered. Instead, a charm parton density in the hadron, c{x, Q 2 ), has to 
be introduced. This density obeys the usual massless NLO DGLAP evolution equation with the boundary condition 
c(x, Q 2 = Qq) = where Qq ~ to 2 . So, we may say that, within the VFNS, the charm density arises perturbatively 
from the g — > cc evolution. 

In the VFNS, the treatment of the charm depends on the values chosen for Q 2 . At low Q 2 < Qq, the production 
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cross sections are described by the light parton contributions (m, d, s and g). The charm production is dominated by 
the GF process and its higher order QCD corrections. At high Q 2 ^> to 2 , the charm is treated in the same way as 
the other light quarks and it is represented by a charm parton density in the hadron, which evolves in Q 2 . In the 
intermediate scale region, Q 2 ~ to 2 , one has to make a smooth connection between the two different prescriptions. 

Strictly speaking, the perturbative charm density is well defined at high Q 2 3> m 2 but does not have a clean 
interpretation at low Q 2 . Since the perturbative IC originates from resummation of the mass logarithms of the 
type a" \n n (Q 2 /m 2 ), it is usually assumed that the corresponding PDF vanishes with these logarithms, i.e. for 
Q 2 < Qq ~ to 2 . On the other hand, the threshold constraint W 2 = (q + p) 2 = Q 2 (l/x — 1) > 4m 2 implies that 
Qo is not a constant but "live" function of x. To avoid this problem, several solutions have been proposed (see e.g. 
Refs. [H, SH|). In this paper, we use the so-called ACOT(x) prescription [43| which guarantees (at least at Q 2 > to 2 ) 
the correct threshold behavior of the heavy-quark-initiated contributions. 

Within the VFNS, the (^-independent charm production cross sections have three pieces: 

a 2 (x, A) = (72,gf{x, A) - o- 2i sub(x, A) + a 2 ,Qs(x, A), (52) 

where the first and third terms on the right hand side describe the usual (unsubtracted) GF and QS contributions 
while the second (subtraction) term renders the total result infra-red safe in the limit to 2 — > 0. The only constraint 
imposed on the subtraction term is to reproduce at high energies the familiar MS partonic cross section: 

Urn [& 3>g (z, A) - a 2 ,su B (z, A)] = cr^(z). (53) 

Evidently, there is some freedom in the choice of finite mass terms of the form X n (with a positive n) in <J2,sub{z, A). 
For this reason, several prescriptions have been proposed to fix the subtraction term. As mentioned above, we use 
the so-called ACOT(x) scheme [tjjj]. 

According to the ACOT(x) prescription, the lowest order (^-independent cross section is 



T 2 (x,X) = I dzg(z,fi F ) 



l 

2 



&Sl(x/z, A)^lnga B (x/z) P^ c (*/*) 



+ <7 B (x)c + (x,/i F ), (54) 



X 



9 



where Pg^lc is the LO gluon-quark splitting function, Pg°l c (() = [(1 — () 2 + C, 2 ] / 2, and the LO GF cross section 
is given by Eqs. (|39jl . Remember also that x = x (l + 4A) and c+(£, ^f) — c{(, Mf) + c(£, /xp)- 

The asymptotic behavior of the subtraction terms is fixed by the parton level factorization theorem. This theorem 
implies that the partonic cross sections d<7 can be factorized into process-dependent infra-red safe hard scattering 
cross sections da, which are finite in the limit to — > 0, and universal (process-independent) partonic PDFs / _ >j and 
fragmentation functions d n ^Q: 

d<5-(7* +a^Q + X)=Y^ fa^iiO ® da( 7 * + i -> n + X) ® d n ^ Q (z). (55) 

In Eq. (|55p . the symbol ® denotes the usual convolution integral, the indices a,i,n and Q denote partons, pi = Qp a 
and pq = zp n . All the logarithms of the heavy-quark mass (i.e., the singularities in the limit to — > 0) are contained 
in the PDFs f a —n and fragmentation functions d n -»Q while da are IR-safe (i.e., are free of the In to 2 terms). The 
expansion of Eq. (|55p can be used to determine order by order the subtraction terms. In particular, for the LO GF 
contribution to the charm leptoproduction one finds Q 

°k?suB (*. In (A4/™ 2 )) = fg^c (C, In {^ F /m 2 )) ® a^ QS (z/0, (* - 2, L, A, I), (56) 

where /gH c (C; m (/^f/™ 2 )) — ( a s/27r) In (/x|./m 2 ) Pg^l c (C) describes the charm distribution in the gluon within the 
MS factorization scheme. 

One can see from Eq. (|56p that the azimuth-dependent GF cross sections o a,gf and oifiF don't have subtraction 
terms at LO because the lowest order QS contribution is (^-independent. For this reason, the cos 2(p-dependence 
within the VFNS has the same form as in the FFNS: 

l 

af°'(x,\) = / Azg(z 1 LL F )a [ ^ g {xlz, A) . (57) 
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FIG. 8: Azimuthal asymmetry parameter A(x, A) in the VFNS at several values of A. The following contributions are plotted: 
GF (LO) (solid curves), GF (LO) -SUB (LO) +QS (LO) with the CTEQ5M set of PDFs (dotted curves) and GF (LO) -SUB (LO) + 
QS (LO) with the CTEQ4M set of PDFs (dashed curves). 



Fig.[8]shows the ACOT(x) predictions for the asymmetry parameter A(x, A) at several values of variable A: A -1 = 
1, 4, 10, 20, 50 and 100. For comparison, we plot also the LO GF predictions (solid curves). In the ACOT(x) case, we 
consider the CTEQ5M (dotted lines) and CTEQ4M (dashed curves) parametrizations of the gluon and charm densities 
in the proton. Corresponding values of the charm quark mass are m c — 1.3 GeV [40J (for the CTEQ5M PDFs) and 
m c = 1.6 GeV [44] (for the CTEQ4M PDFs). The default value of the factorization scale is Hf = V™ 2 + Q 2 - 

One can see from Fig. |S] the following properties of the azimuthal asymmetry, A(x, A), within the VFNS. Contrary 
to the nonperturbative IC component, the perturbative one is significant practically at all values of Bjorken x and 
Q 2 > to 2 . The perturbative charm contribution leads to a sizeable decreasing of the GF predictions for the cos2iy9- 
asymmetry. In the ACOT(x) scheme, the IC contribution reduces the GF results for A(x, A) by about 30%. The 
origin of this reduction is straightforward: the QS component is practically cos 295-independent. 
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The ACOT(x) predictions for the asymmetry depend weakly on the parton distribution functions we use. It is seen 
from Fig. E] that the CTEQ5M and CTEQ4M sets of PDFs lead to very similar results for A(x, A). Note that we give 
the CTEQ5M predictions at low x only because of irregularities in the CTEQ5M charm density at large x. 

We have also analyzed how the VFNS predictions depend on the choice of subtraction prescription. In particular, 
the schemes proposed in Refs. [H, EE) have been considered. We find that, sufficiently above the production threshold, 
these subtraction prescriptions reduce the GF results for the asymmetry by approximately 30 ~ 50%. 

One can conclude that impact of the perturbative IC on the cos 2ip asymmetry is essential in the whole region of 
Bjorken x and therefore can be tested experimentally. 



In the present paper, we consider the azimuthal dependence in charm leptoproduction as a probe of the IC content 
of the proton. Our analysis is based on the fact that the GF and QS components have strongly different cos2<^- 
distributions. This fact follows from the NLO calculations of both parton level contributions. In the framework of 
the FFNS, we justify the most remarkable property of the hadron level azimuthal cos2</? asymmetry: the combined 
GF+QS predictions for A(x,Q 2 ) are perturbatively and parametrically stable. The nonperturbative IC contribution 
(resulting from the five-quark \uudcc) component of the proton wave function) is practically invisible at low a:, but 
affects essentially the GF predictions for the asymmetry at large Bjorken x. We conclude that measurements of the 
cos 2tp asymmetry at large x could directly probe the nonperturbative intrinsic charm. 

Within the VFNS, charm density originates perturbatively from the g — > cc process and obeys the DGLAP evolution 
equation. Presently, charm densities are included practically in all the global sets of PDFs like CTEQ and MRST. 
Our analysis shows that these charm distribution functions reduce dramatically (by about 1/3) the GF predictions 
for A(x, Q 2 ) practically at all values of x. For this reason, the perturbative IC contribution can easily be measured 
using the azimuthal cos 2c/?-distributions in charm leptoproduction. 

The VFN schemes have been proposed to resum the mass logarithms of the form a™ ln n (Q 2 /m 2 ) which dominate 
the production cross sections at high energies. Evidently, were the calculation done to all orders of a s , the VFNS and 
FFNS (without nonperturbative IC) would be exactly equivalent. There is a point of view advocated in Refs. @, @] 
that, at high energies, the perturbative series converge better within the VFNS than in the FFNS. There is also 
another opinion [46l . |47| that the above logarithms do not vitiate the convergence of the perturbation expansion so 



that a resummation is, in principle, not necessary. Our analysis of the azimuthal dependence in leptoproduction 
indicates an experimental way to resolve this problem. First, contrary to the production cross sections, the azimuthal 
cos 2iy9-asymmetry is well defined numerically in pQCD. Second, sufficiently above the production threshold (i.e., at 
small enough Bjorken x), the LO VFNS predictions for the cos 2^-asymmetry differ by more than 30% from the 
corresponding FFNS ones. Third, nonperturbative contributions (like the intrinsic gluon motion in the target) can't 
compensate for this difference at non-small Q 2 where the VFNS is expected to be adequate. Therefore measurements 
of the azimuthal distributions in charm leptoproduction would make it possible to clarify the question whether the 
VFNS perturbative series for A(x 1 Q 2 ) converges better than the FFNS one. 



We thank S.J. Brodsky for stimulating discussions and useful suggestions. We also would like to acknowledge 
interesting correspondence with I. Schienbein. This work was supported in part by the ANSEF grant 04-PS-hcpth- 
813-98 and NFSAT grant GRSP-16/06. 

APPENDIX A: VIRTUAL AND SOFT CONTRIBUTIONS TO THE QUARK SCATTERING 

In this Appendix we reproduce some results of Hoffmann and Moore for the ^-independent QS cross sections, 
and correct two misprints uncovered in Rcf. [25| . We work in four dimensions, in the Fcynman gauge and use the 
on-mass-shell renormalization scheme. We compute the absorptive part of the Feynman diagram (which is free of the 
UV divergences) and then restore the real part using the appropriate dispersion relations. 

In the on-mass-shell scheme, the renormalized fermion self-energy vanishes like (pq — m) 2 which means that the 
second and third diagrams in Fig. do not contribute to the cross section when the external quark legs are on-shell, 
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Pq — > m. The first graph in Fig.[2J; describes the NLO corrections to the quark-photon vertex function: 



g(Q 2 

2m 



(Al) 



where = | (7^7^ — 7^7^) while / (Q 2 ) and g (Q 2 ) are the quark electromagnetic formfactors. At the lowest order 



The virtual lepton-quark cross section, ct^, is obtained from the interference term between the virtual and the 
Born amplitude. The result can be written in terms of the electromagnetic formfactors as: 



d 2 ~ v 



IQ 



dzAQ 2 



zQ 2 



S(i-z){[i + (i-y) 2 - 2Az V] / (Q 2 ) + y 2 g (Q 2 ) } 



Taking into account the definition of the HM cross sections, and a^ L \ given by Eqs. 
that corresponding virtual parts are: 



a (2 v (z,Q 2 ) = 2S(l-z)f^ ( Q 2) 



a[ L v \z,Q 2 ) = -5(l-z)g^ (Q 2 ) 



(A2) 

and (fT5|) . we find 
(A3) 



where /« (Q 2 ) and ff W (Q 2 ) are the NLO corrections to the electromagnetic formfactors. For the NLO HM cross 

sections, crjy and u[y , we use exactly the same notations as in Ref. (25j |. 

In the on-mass-shell renormalization scheme, the renormalized vertex correction vanishes as the photon virtuality 
goes to zero, /W (0) = 0. This is a consequence of the Ward identity and the fact that the real photon field (like 
the massive fermion one) is unrenormalized in first order QCD. To satisfy the condition /W (0) = automatically, 
we should use for /W (g 2 ) the dispersion relation with one subtraction. The second formfactor, g^ (q 2 ), has no 
singularities. For this reason, we use for g^ (g 2 ) the dispersion relation without subtractions: 



(. 



n = q - 



dt Im fW(t) 



t(t- 



■ ioy 



,(i) 



00 



dt Im ffWft) 
t - q 2 - iO 



(A4) 



Calculating the imaginary parts of the formafactors and restoring their real parts with the help of Eqs. (|A4[) yields 



/w (g 2 ) = 



,(i) 



(Q 2 ) = 



1 + 2A 
VI + 4A 



In 1 



1 + 2A 
vT + 4A 



a s A In r 
— Of ■ — ^=. 
7T VI +4A 



12 



ln^-1 



Lia(-r) + — + 



■ In 2 r 



— In 7' In — — 

2 A 



1 lnr 



4 V1 + 4A 



(A5) 
(A6) 



In Eqs. |[A"5]) and (|A"6|l . G> = (iV| - l)/(2N c ), where iV c is number of colors, and r is defined by Eq. (27)). Taking 
into account that A^ 1 ^ = (a s /a e m ) Cf Aq^, we see that Eqs. (IA5[ IA6I) reproduce the textbook QED results. 

It is now straightforward to obtain the virtual contribution to the longitudinal cross section. Combining Eqs. 
and (|A6|) yields: 



4v\z,Q 2 ) 



a s A lnr 

— Gfo(1 — z) ; =. 



(A7) 



Comparing the above result with the corresponding one given by Eq. (39) in Ref. [251 ], we see that the HM expression 
for a^y has opposite sign. Note also that this typo propagates into the final result for cr^ given by Eq. (52) (25| . 

Calculation of the bremsstrahlung contribution to the longitudinal cross section, cr^g(z, Q 2 ), is also straightforward. 
We coincide with the HM result for a[^(z, Q 2 ) given by Eq. (49) in Ref. [25| . However there is one more misprint in 
the HM expression for (j[ L ^: the r.h.s of Eq. (52) (2f| should be multiplied by z. 

In the case of cr^(z,Q 2 ), the situation is slightly more complicated due to the need to take into account the IR 
singularities. One can see from Eq. (|A5| that /W (Q 2 ) has an IR divergence which is regularized with the help of an 
infinitesimal gluon mass m g . This singularity is cancelled when one adds the so-called soft contribution originating 
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from the real gluon emission. For this purpose we introduce another infinitesimal parameter Sz, (m g /m) <C Sz <C 1. 

(2) 

The full brcmsstrahlung contribution, cr^, can then be splitted into the soft and hard pieces as follows: 



a (2 J oft (z, Q 2 ) = 6(z + Sz- l)a^(z, Q 2 ), 



<r\Ld& Q 2 ) = 0(l-z- 6z)a[ 2 >(z, Q 2 ), 



(A8) 



where 6(1 — z — Sz) is the Heaviside step function. The soft cross section should be calculated in the eikonal 
approximation, p g — > 0, taking into account the infinitesimal gluon mass m g . As a result, the sum of the virtual 
and soft contributions is IR finite: 



r (2) 



r (2) 
'lsoft 



— C F 5{l-z){ -2\n{5z 



1 + 2A 
VI + 4A 



1 + 2A 
V1 + 4A 



: lnr 



2 In A - 1 



V1 + 4A 



lnr 



Li 2 (r 2 ) + 2Li 2 (-r) + - In 2 r - 2 In r - In r In A + 2 In r ln(l + 4A) 



(A9) 



(2) 

Adding to the above expression the hard cross section c\i! ard defined by Eq. 
the full result for er^ given by Eq. (51) in Ref. (2o| . 



we reproduce in the limit Sz — > 



APPENDIX B: NLO SOFT-GLUON CORRECTIONS TO THE PHOTON-GLUON FUSION 

This Appendix provides an overview of the NLO soft-gluon approximation for the photon-gluon fusion mechanism. 
We present the final results for the parton level cross sections to the next-to-leading logarithmic (NLL) accuracy. 
More details can be found in Refs. [Hj (H [H . 

To take into account the NLO contributions to the GF mechanism, one needs to calculate the virtual 0(a em a 2 ) 
corrections to the Born process and the real gluon emission: 

7 *(g) + g(k g ) -> Q( PQ ) + Q(p-)+g( Pg ). (Bl) 
The partonic invariants describing the single-particle inclusive (1PI) kinematics are 

s' = 2q-k g = s + Q 2 = £S', h = {k g - p Q f - m 2 = £T U 

S4 = s' + ti+Ui, ui = (q -Pq) 2 - m 2 = Ui, (B2) 

where £ is defined by k g = £p and S4 measures the inelasticity of the reaction (|B1|) . The corresponding 1PI hadron 
level variables describing the reaction ([T]) are 

5' = 2q-p = S + Q 2 1 T^ip-pQ) 2 -m 2 , 

S 4 = S' + T 1 + U 1 , U 1 = (q- p Q f - m 2 . (B3) 

The exact NLO calculations of the unpolarized heavy quark production in 73 [48|,[4i|, 7*3 (38|, and gg [50L [5ll[52l. l53j 
collisions show that, near the partonic threshold, a strong logarithmic enhancement of the cross sections takes place in 
the collinear, p g T — * 0, and soft, p g — > 0, limits. This threshold (or soft-gluon) enhancement has universal nature in 
the perturbation theory and originates from incomplete cancellation of the soft and collinear singularities between the 
loop and the brcmsstrahlung contributions. Large leading and next-to-leading threshold logarithms can be resummed 
to all orders of perturbative expansion using the appropriate evolution equations [54l [55l[56l|. The analytic results 
for the resummed cross sections are ill-defined due to the Landau pole in the coupling strength a s . However, if one 
considers the obtained expressions as generating functionals of the perturbative theory and re-expands them at fixed 
order in a s , no divergences associated with the Landau pole are encountered. 

Soft-gluon resummation for the photon-gluon fusion has been performed in Ref. and checked in Refs. [13, HH . 
To NLL accuracy, the perturbative expansion for the partonic cross sections, d 2 <7fc iff /diidwi (k — T, L, A, I), can be 
written in a factorized form as 
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with the Born level distributions B®° rn given by 
B*°™{s'M,u x ) 
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(B5) 



1 - 



2Q 2 



2m 2 s' 



Note that the functions if( n )(s' , ti,ux) in Eq. (|B4|) originate from the collinear and soft limits. Radiation of soft 
and collinear gluons does not affect the transverse momentum of detected particles and therefore the azimuthal angle 
tp. For this reason, the functions K^ n \s' ,ti,u\) are the same for all helicity cross sections &k, g (k — T, L, A, I). At 
NLO, the soft-gluon corrections to NLL accuracy in the MS scheme are 
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(B6) 



where we use /i = fiF = (Ir- In Eq. (|B6|) . Ca = N c and Cf = (A 2 — 1)/(2N C ), where N c is number of colors, while 
Lp = (1 — 2m 2 /s){ln[(l — /3)/(l + /3)]+i7r} with j3 = \J\ — Am 2 /s. The single-particle inclusive "plus" distributions 
are defined by 
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(B7) 



For any sufhciently regular test function ft.(s 4 ), Eq. (IB7|) gives 
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y ds 4 /l(s 4 ) 
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(B8) 



In Eq. (|B6[) . we have preserved the NLL terms for the scale-dependent logarithms too. Note also that the results 
(IB5|) and (|B6p agree to NLL accuracy with the exact 0(a em a 2 ) calculations of the photon-gluon cross sections ar. g 
and &L, g given in Ref. [38j ] . 

To investigate the scale dependence of the results (|B4HB6J) . it is convenient to introduce for the fully inclusive 



(integrated over t\ and u\) cross sections, &k. g (k = T, L, A, I), the dimensionless coefficient functions cj^ defined 

by Eq. Concerning the NLO scale- independent coefficient functions, only d^'g and c^ 1 ^ are known exactly 

[38l[57j. As to the u-dependent coefficients, they can by calculated explicitly using the evolution equation: 

l 

dCa k , g (z/(,Q 2 ,n 2 )P gg (C), (B9) 



d<7 fe , g (z,Q 2 ,/r) 



d In /i 2 



Cmin 



where z = Q / s 1 , Q 



z(l + 4A), ak, g (z, Q 2 ,fj) are the cross sections resummed to all orders in a s and P gg (C) is 



the corresponding (resummed) Altarelli-Parisi gluon-gluon splitting function. Expanding Eq. (|B9j) in a s , one can find 

MM 



4 1 3 1) (^A) = i ^ / 6£[h6{l-Q- P£\Q]$f(z/C,\), 



(BIO) 
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where 62 = (HCU — 2ny)/12 is the first coefficient of the /3(a s )-function expansion and n/ is the number of active 
quark flavors. The one-loop gluon splitting function is: 

PW(0 = lim { + ^ + C(l - 0) e(l-C-e) + 6(1 - lnej + fo 2 5(l - 0- (Bll) 

With Eq. (|B10I) in hand, it is possible to check the quality of the NLL approximation against exact answers. As 
shown in Ref. [2l|, the soft-gluon corrections reproduce satisfactorily the threshold behavior of the available exact 
results for A ^ 1. Since the gluon distribution function supports just the threshold region, the soft-gluon contribution 
dominates the photon-hadron cross sections Ok,GF (k = T, L, A, I) at energies not so far from the production threshold 
and at relatively low virtuality Q 2 < m 2 . 



APPENDIX C: NONPERTURBATIVE IC AND RELEVANT EXPERIMENTAL FACTS 

The most clean probe of the charm quark distribution function (both perturbative and nonperturbative) is the 
semi-inclusive deep inelastic lepton-proton scattering, Ip — > VcX. To measure the nonperturbative IC contribution, 
one needs data on the charm production at sufficiently large Bjorken x. The only experiment which has investigated 
the large x domain is the European Muon Collaboration (EMC) [H^] where the decay lepton spectra have been used to 
detect the produced charmed particles. In Ref. [H^], a re-analysis of the EMC data on F£(x, Q 2 ) have been performed 
using the NLO results for both GF and QS components. The analysis [59| shows that a nonperturbative intrinsic 
charm contribution to the proton wave function of the order of 1% is needed to fit the EMC data in the large x region. 
This value of the nonperturbative IC is consistent with the estimates based on the operator product expansion [{|. 
Note however that the EMC data are of limited statistics and, for this reason, more accurate measurements of charm 
leptoproduction at large x are necessary. 

It is also possible to extract useful information on the IC from diffractive dissociation processes such asp-> p J/ i/j 
on a nuclear target. Comprehensive measurements of the pA — > J/tpX and irA — ► J/tpX cross sections have been 
performed in the fixed target experiments NA3 at CERN and E886 at FNAL [6l| . According to the arguments 
presented in Refs. [H, HHHi]], the IC contribution is predicted to be strongly shadowed in the above reactions that 
is in a complete agreement with the observed nuclear dependence of the high Feynman xf component of the J/ijj 
hadroproduction. 

A non- vanishing five-quark Fock component \uudcc) leads to the production of open charm states such as h c (cud) 
and D~ (cd) with large Feynman xp. This may occur either through a coalescence of the valence and charm quarks 
which are moving with the same rapidity or via hadronization of the produced c and c. As shown in Refs. [rja . [6^ |. 
a model based on the nonperturbative intrinsic charm naturally explains the leadi ng p article effect in the pp — ► DX 
and pp — > A c X processes that has been observed at the ISR [67] and Fermilab [68|, l69( . 

As to the high-Xf hadroproduction of open bottom states like Ab(bud), corresponding cross sections are predicted 
to be suppressed as m 2 /m 2 ~ 1/10 in comparison with the case of charm production. Evidence for the forward Af, 
production in the pp collisions at the ISR energy was reported in Refs. [70L l71|. 

Rare seven-quark fluctuations of the type \uudcccc) in the proton wave function can lead to the production of two 
J/tjj [lH or a double-charm baryon state at large Xp and low pt- Double J/ip events with a high combined xp > 0.5 
have been detected in the NA3 experiment An observation of the double-charmed baryon 2+ (3520) with mean 
(x F ) ~ 0.33 has been reported by the SELEX collaboration at FNAL [zU. 
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